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Section  I 


INTRODUCTION 


Air  Force  analysis  of  operational  experience  with  aircraft 
fleets  over  the  past  several  years  has  led  to  the  establishment 
of  damage  tolerance  criteria  for  metal  airframes  [1]  which  require 
the  designer  to  verify  the  structural  integrity  of  flight-safety- 
critical  details  based  on  fracture  mechanics  calculations.  These 
calculations  are  distinct  from  and  serve  a different  purpose  than 
fatigue  calculations  based  on  Miner’s  rule  or  similar  damage  accumu- 
lation hypotheses.  The  damage  tolerance  philosophy  assumes  the 
presence  of  small  cracks  in  the  airframe  at  the  time  the  airplane 
is  brought  into  service,  and  seeks  to  protect  the  structure  against 
failures  which  might  result  from  the  actual  presence  of  these  cracks. 
The  methodology  is  based  upon  integration  of  an  empirical  crack  growth 
rate  relationship  over  the  aircraft  life  history,  as  represented  by 
loads  in  the  vicinity  of  the  safety-critical  parts,  and  upon  the 
establishment  of  a fracture-critical  crack  size  for  each  part 
corresponding  to  a maximum  expected  load. 

Both  the  crack  growth  calculations  and  the  determination  of 
fracture-critical  crack  sizes  require  the  application  of  linear 
elastic  fracture  mechanics  to  compute  stress  intensity  factors  at 
the  tip  of  the  assumed  sharp  crack.  Analytical  solutions  for  stress 
intensity  factors  can  be  obtained  for  crack  tips  remote  from  other 
geometrical  features  of  the  structure  boundary.  These  solutions 
may  be  extended  to  more  complicated  geometries  by  the  techniques 
of  conformal  mapping  and  boundary  collocation  [2,3].  However,  the 
complexities  found  in  aircraft  structural  details  usually  involve 
the  proximity  of  two  different  types  of  boundary  geometry  (e.g.,  a 
circular  fastener  hole  near  a linear  edge)  as  well  as  load  distribu- 
tions which  are  more  complicated  than  uniform  stress,  linearly  vary- 
ing stress,  or  a point  load.  Thus,  stress  intensity  factor  solutions 
by  the  traditional  methods  are  often  inconvenient,  and  may  be 
prohibitive  in  some  cases. 
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The  above  problems  are  also  encountered  in  conventional  stress 
analysis.  The  finite  element  methods  have  assumed  a dominant  role 
in  conventional  stress  analysis  over  the  past  decade  because  of 
their  inherent  ability  to  accommodate  complex  boundaries  and  load 
distributions  in  a highly  organized  fashion  which  is  also  convenient 
for  the  analyst.  For  this  same  reason,  the  finite  element  methods 
are  now  beginning  to  be  applied  extensively  to  the  field  of  fracture 
mechanics  analysis,  A recent  paper  by  Pian  [4]  reviews  several 
alternate  approaches  to  finite-element  fracture  mechanics  and  com- 
pares their  relative  merits.  The  work  reviewed  in  the  present 
report  has  been  based  on  the  assumed-stress  hybrid  approach  [5] . 

The  key  feature  in  applying  the  hybrid  finite  element  method 
to  fracture  mechanics  is  the  special-purpose  crack-tip  element 
shown  in  Fig.  1.  This  element  possesses  an  assumed  displacement 
distribution  along  its  outer  boundary  which  is  compatible  with  the 
boundary  displacement  assumption  for  conventional  quadrilateral 
elements.  The  stress  assumption,  which  characterizes  behavior  within 
the  element,  includes  the  crack-tip  singularity  and  several  additional 
nonsingular  terms,  all  of  which  are  obtained  from  a complex  variable 
elasticity  solution  for  the  material  region  near  the  crack  tip  [6] . 

In  effect,  what  has  been  done  is  to  combine  the  classical  elasticity 
solution  for  a crack-tip  stress  singularity,  a simple  solution  in 
the  absence  of  neighboring  boundaries  and  applied  loads,  with  the 
flexibility  of  transitioning  to  these  boundaries  and  loads  in  the 
piece-wise  fashion  permitted  by  the  finite  element  methods.  From 
a practical  standpoint,  the  hybrid  approach  also  proves  to  be  con- 
venient because  it  permits  the  creation  and  solution  of  finite  ele- 
ment models  using  the  well-known  Matrix  Displacement  Method.  The 
special  crack-tip  element  is  in  this  respect  just  another  element 
with  nodes  to  be  coupled  to  the  model.  Its  only  unusual  features 
are  its  shape,  its  possession  of  nine  instead  of  four  nodes,  and 
the  fact  that  the  matrix  "stress"  analysis  procedure  one  normally 
associates  with  a conventional  element  leads  in  this  case  to  stress 
intensity  factors  instead  of  stresses. 


2 


Naturally,  there  is  a price  to  pay  for  the  convenience  of  the 
hybrid  method.  Crack  tips  cannot  be  located  near  other  geometrical 
details  with  complete  disregard  for  the  crack-tip  position  relative 
to  the  boundaries  of  the  special  element.  One  is  thus  forced  to 
accept  a compromise  between  solution  accuracy  and  the  cost  of  the 
analysis,  parameters  which  increase  as  the  finite  element  model 
becomes  more  detailed.  The  present  work  has  focussed  on  analysis 
of  several  typical  aircraft  structural  details  and  has  shown  that 
the  hybrid  approach  permits  a reasonable  compromise  between  cost 
and  accuracy.  Subsequent  sections  of  this  report  review  the  scope 
of  the  program,  tests  for  accuracy,  computation  cost  experience, 
and  the  detailed  results  obtained  from  parametric  analyses  of  some 
of  the  structural  details. 
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Section  II 


PROGRAM  OVERVIEW 

The  major  objectives  of  the  program  were  to  provide  numerical 
analysis  codes  to  compute  the  mode  I and  mode  II  stress  intensity 
factors  K^- , for  four  general  classes  of  structural  details: 

an  attachment  lug  with  a bearing  hole,  a rectangular  stiffened  or 
unstiffened  panel  with  an  open  hole  which  could  be  centered  or  off- 
set, a single  row  of  fastener  holes  near  one  edge  of  a panel,  and 
a double  staggered  row  of  fastener  holes  near  one  edge  of  a panel. 

In  all  cases,  cracks  were  to  be  assumed  emanating  radially  from  one 
of  the  holes:  a single  crack  for  the  attachment  lug  and  either  one 

crack  or  two  diametrically  opposed  cracks  in  the  other  configurations. 
The  angular  locations  of  the  cracks  around  the  hole  were  to  be  vari- 
able, and  selection  of  the  damaged  hole  was  to  be  a parameter  in  the 
cases  of  the  multi-hole  details.  The  general  configurations  are 
illustrated  schematically  in  Fig.  2. 

The  above  objectives  were  achieved  by  programming  individual 
parametric  analysis  codes  for  each  configuration,  using  the  ASRL 
FEABL-2  software  [7]  as  a base.  Three  interim  technical  reports 
[8,9,10]  discuss  in  detail  the  construction,  verification  and 
initial  demonstration  of  each  computer  program.  Additional  infor- 
mation relating  to  accuracy  analysis  and  the  general  capabilities 
of  hybrid  crack  elements  has  also  been  published  separately  [11,12, 
13,14].  An  additional  objective  to  study  the  effect  on  and  K 
of  an  interference-fit  fastener  in  the  rectangular  panel  (Fig.  2B) 
could  not  be  attempted  within  the  scope  of  the  program.  The  load- 
ing applied  to  the  fastener  hole  row  details  (Fig.  2C  and  2D)  was 
changed  from  uniform  tension  on  the  neighboring  edge  to  individual 
bearing  loads  on  each  hole  during  the  program,  when  the  latter 
situation  proved  to  be  more  interesting  in  the  light  of  results 
obtained  from  analyses  of  the  rectangular  panel  configuration. 

It  should  be  observed  that  the  stress  intensity  factors  K_[ , 

K j will  vary  with  the  material  Poisson  ratio  v and  will  have 
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different  values  in  plane  stress  and  plane  strain  when  the  structure 
incorporates  a crack  tip  near  a constrained  boundary  or  when  it  is 
otherwise  multiply  connected,  as  in  the  present  cases.  It  has  been 
remarked  [3]  that  the  thickness  of  a structure  relative  to  its 
external  dimensions  determines  whether  plane  stress  or  plane  strain 
analysis  is  to  be  conducted,  irrespective  of  whether  plane  strain 
conditions  are  achieved  locally  at  the  crack  front.  The  plane  stress 
state  was  considered  to  be  of  primary  interest  in  the  present  work, 
and  the  computer  codes  have  been  programmed  accordingly.  In  all  of 
the  analyses  discussed  subsequently,  results  were  obtained  based  on 
the  general  properties  for  aluminum  alloys: 

E = 107  psi  v = 0 . 3 

The  elastic  properties  E,  v are  input  parameters,  but  the  sensitivity 
of  K , solutions  to  v was  not  studied  in  the  present  program. 

Where  possible,  the  finite  element  solutions  for  stress  intensity 
factors  were  compared  with  independently  obtained  analytical  or  bound- 
ary collocation  solutions  [3].  However,  since  independent  solutions 
are  not  available  for  most  of  the  configurations  and  crack  locations 
considered  in  the  present  work,  accuracy  assessments  were  also  made 
based  upon  the  performance  of  individual  types  of  finite  elements 
under  varying  conditions.  The  isoparametric  4-node  quadrilateral 
was  used  as  the  basic  building-block  in  the  various  finite  element 
models.  Its  performance  has  been  well  characterized,  and  it  is 
known  to  give  accurate  displacement  solutions  in  the  presence  of 
gradients  if  the  element  shape  is  kept  close  to  a square  [15] . 

The  mesh-generation  portions  of  the  computer  codes  were  programmed 
accordingly,  and  attention  was  focussed  upon  the  behavior  of  the 
hybrid  crack  element  (PCRK59).  The  element  was  found  to  have  a 
somewhat  better  performance  with  respect  to  shape  distortions, 
using  accuracy  of  the  computed  value  as  a measure.  The  key  error 
parameter  was  found  to  be  the  position  of  the  crack  tip  within  the 
element  [8,13].  The  latter  effect  is  illustrated  in  Fig.  3,  which 
plots  the  error  in  computed  stress  intensity  for  the  problem  of  a 
finite-width  strip  with  symmetrical  edge  cracks  loaded  in  uniform 
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tension.  As  can  be  seen  in  the  figure,  serious  errors  are  introduced 
if  the  crack  tip  approaches  the  opposite  edge  of  the  PCRK59  element. 
This  effect  is  attributed  to  the  tendency  of  severe  local  gradients 
to  "spill  over"  into  the  neighboring  quadrilateral  elements,  which 
are  not  equipped-  to  handle  such  situations.  The  results  also  seem 
to  indicate  that  reasonably  accurate  solutions  are  obtained  when 
the  crack  tip  is  near  a lateral  boundary  of  the  structure  (also 
the  PCRK59  boundary),  i.e.,  for  a/b  = 0.2  in  Fig.  3.  However, 
additional  results  discussed  subsequently  cast  some  doubt  on  this 
finding. 

The  attachment  lug  analysis  was  programmed  with  the  quadrilateral 
and  PCRK59  elements  as  the  basic  building  blocks.  However,  initial 
attempts  to  design  an  adequate  mesh  for  the  rectangular  panel  led  to 
a mesh-grading  problem  which  was  finally  solved  by  creating  another 
building  block.  As  shown  in  Fig.  4A  and  4B,  the  amount  of  refinement 
required  near  the  fastener  hole  demands  either  a wasteful  use  of 
nodes  and  elements  throughout  the  panel  or  a rapid  grading  to  a very 
coarse  mesh  remote  from  the  hole.  While  the  first  course  of  action 
would  have  resulted  in  a costly  analysis,  the  second  would  have  led 
to  inaccuracy  in  the  displacement  solution  because  the  mesh  would 
have  contained  constant-strain  triangle  elements.  Instead,  a third 
alternative  was  chosen  by  using  the  hybrid  method  to  create  a 
special-purpose  "hole  element"  which  occupies  the  region  between 
the  square  and  outer  circular  boundaries  in  Fig.  4B.  Since  the 
essence  of  the  hybrid  approach  to  special-purpose  elements  is  to 
mimic  known  behavior,  mid-edge  nodes  were  included  on  the  hole  ele- 
ment outer  boundary  to  allow  in  an  average  sense  for  the  type  of 
displacement  field  expected  in  the  vicinity  of  an  open  hole,  and 
the  assumed  stress  field  was  adapted  from  a classical  elasticity 
solution  for  this  type  of  region. 

Extensive  tests  of  the  hole  element  with  and  without  inner 
circular  rings  of  quadrilaterals  verified  the  element* s ability  to 
reproduce  accurate  stress  fields  for  a hole  unloaded,  subjected  to 
uniform  pressure  (interference-fit)  or  subjected  to  a bearing  load; 
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tests  with  PCRK59  elements  in  place  in  the  inner  rings  also  verified 
the  ability  of  the  finite  element  model  to  compute  stress  intensity 
factors  [9].  Figure  4C  illustrates  a typical  solution  obtained  from 
the  rectangular  panel  program,  with  the  fastener  hole  on  center  and 
without  edge  stiffeners.  This  "butterfly”  plot  summarizes  K^_, 
as  functions  of  the  polar  angle  to  the  crack.  The  solution  for 
was  verified  for  single  cracks  at  0°  and  pairs  of  cracks  at  0°,  180° 
by  comparison  with  independent  solutions  [3] . 

The  results  reported  previously  for  panel  program  performance 
[9]  were  considered  to  be  incomplete,  partly  because  the  independent 
solution  which  was  then  available  did  not  account  for  finite-width 
effects,  and  partly  because  sensitivity  to  the  shape  of  the  hole 
element  had  not  been  thoroughly  investigated.  The  study  of  both 
of  these  performance  aspects  has  since  been  completed  after  receipt 
of  some  revisions  to  Ref.  3 which  included  finite-width  solutions 
for  cracks  at  0°  and  180°  at  a hole  in  a tension  panel.  Figure  5 
illustrates  the  effect  of  the  ratio  of  crack  size  to  hole  radius 
a/R  on  solution  accuracy.  The  error  in  is  seen  to  range  between 
-9  to  +3  percent,  with  the  best  performance  occuring  for  0.05  _<  a/R 
£ 1.2.  The  upper  limit  is  sufficient  to  extend  the  range  of  good 
solutions  to  regions  in  which  the  fastener  hole  begins  to  lose  its 
influence  on  the  crack  tip.  For  very  small  crack  sizes  (a  = 0.005 
inch,  a/R  - .0158)  the  solution  error  is  observed  to  be  significantly 
greater  than  was  previously  estimated  (Fig.  3 and  discussion  on  pp.  5-6) 
for  small  cracks.  The  non-monotonic  error  variation  shown  in  Fig.  5 
is  attributed  to  the  manner  in  which  the  mesh  generation  is  programmed 
for  the  interior  rings  of  quadrilaterals  (see  Fig.  4B)  which  surround 
the  fastener  hole.  The  PCRK59  element  is  placed  in  the  first  and 
second,  second  and  third,  or  third  and  fourth  rings  according  to 
the  crack  size,  so  that  the  best  chance  of  placing  the  crack  tip 
at  the  element  center  is  obtained.  This  leads  to  a non-monotonic 
variation  of  the  crack-tip  location  in  the  PCRK59  element,  which 
causes  the  error  behavior. 

The  "shape"  of  the  surrounding  hole  element,  as  expressed  by 
the  ratio  of  its  edge  dimension  to  its  inner  diameter  W/D  , is 
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another  important  performance  parameter.  Results  obtained  from 

earlier  tests  [9]  indicated  that  reasonably  accurate  stresses  and 

stress  intensity  factors  could  be  obtained  over  the  range  1.59  <_ 

W/Dq  _<  6.35.  It  is  desirable  to  extend  this  range  if  possible, 

particularly  the  lower  limit,  in  order  to  permit  the  analysis  of 

other  details  with  closely  spaced  fastener  holes.  Since  the  inner- 

ring  mesh  is  programmed  to  maintain  a fixed  ratio  between  the  fastener 

hole  diameter  D and  the  hole  element  diameter  D : 

o 

D /D  - 2. 52 
o 

the  realizable  center-center  fastener  spacing  is  actually  given  by: 

C/D  =2.52  W/Dq 

Therefore,  some  additional  tests  were  run  to  determine  the  extent 
to  which  the  shape  of  the  hole  element  could  be  distorted.  The 
results,  illustrated  in  Fig.  6,  show  that  surprisingly  accurate 
answers  are  obtained  over  a much  larger  range  than  was  previously 
thought  possible.  The  worst  behavior  (data  points  at  W/Dq  = 1.04, 

1.1)  is  seen  to  be  about  10  percent  error  for  K . Based  on  these 
results,  the  hole  element  is  considered  to  be  useful  over  the  range: 

1.2  < W/D  < 10 
— o — 

The  study  of  sensitivity  to  a/R  (Fig.  5)  was  run  at  W/Dq  - 1.27. 

In  the  limited  parameter  study  which  was  run  previously  to  test 
the  various  panel  program  options  [9],  it  was  found  that  offsetting 
the  hole  from  the  panel  centerline  and/or  placing  lateral  stiffeners 
along  one  or  both  edges  had  very  little  effect  on  the  stress  intensity 
factor  solutions.  This  result  should  be  expected  for  panels  loaded 
by  uniform  tension,  unless  the  crack  tip  is  allowed  to  approach  very 
closely  to  a free  edge  or  a lateral  stiffener.  In  fact,  the  analytical 
solutions  obtained  by  Isida  [16]  and  reproduced  by  finite  element 
analysis  related  to  the  present  work  [14]  indicate  that  the  distance 
between  a crack  (without  fastener  hole)  in  a plate  and  an  edge  stiffener 
must  be  less  than  30  percent  of  the  crack  half-length  to  obtain  more 
than  20  percent  change  in  a solution  which  ignores  the  presence 
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of  the  stiffener.  A similarly  close  approach  is  considered  to  be 
necessary  to  experience  the  effect  of  an  edge  or  stiffener  on  the 
tip  of  a crack  emanating  from  a fastener  hole.  However,  the  panel 
program  was  not  designed  to  handle  these  extremely  close  approaches. 
Therefore,  no  further  parameter  studies  with  this  code  were  conducted. 

The  attachment  lug  and  fastener  row  programs  gave  more  interest- 
ing results  and  were  consequently  exercised  in  extended  parameter 
studies.  The  results  of  these  studies  are  discussed  separately  in 
Sections  3 and  4.  Computing  experience  for  all  of  the  programs  is 
summarized  in  Table  1,  which  compares  typical  total  core  storage 
requirements  and  CPU  times  for  one  complete  set  of  K^,  solutions. 

A "complete  set"  of  solutions  refers  to  '24  to  48  pairs  of  , K 
values  corresponding  to  all  angular  positions  of  the  crack  around 
a fastener  or  bearing  hole,  for  one  given  set  of  dimensions  (hole 
radius,  crack  size,  etc.). 
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Section  III 


ATTACHMENT  LUG  ANALYSIS 

The  attachment  lug  program  automatically  generates  and  analyzes 
the  structural  detail  shown  in  Fig.  7.  Although  the  program  is  able 
to  analyze  a two-material  lug,  with  ideal  bonding  assumed  between  the 
bushing  and  the  remainder  of  the  structure,  this  option  has  not  been 
exercised.  The  other  available  options  (bearing  load  applied  in  one 
of  the  four  directions  indicated  in  the  figure  and  representation  of 
the  load  as  either  a cosine  or  uniform  pressure  distribution)  have 
been  exercised  on  analyses  of  single-material  lug  details.  Detailed 
documentation  of  the  program  and  the  accuracy  tests  associated  with 
it  appear  in  the  first  report  in  this  series  [8]. 

Some  of  the  previous  results  are  reviewed  briefly  here  to  give 
a general  idea  of  the  finite  element  model  and  its  capabilities. 
Figure  8 illustrates  some  of  the  mesh  details  and  an  example  solu- 
tion for  the  case  of  cosine-distributed  tension  bearing  (load  posi- 
tion 1 in  Fig.  7) . The  number  of  quadrilateral  element  divisions 
around  the  bearing  hole  is  an  input  parameter  and  is  usually  speci- 
fied as  24,  32  or  48.  Maximum  shape  distortion  of  the  PCRK59  ele- 
ment occurs  when  24  divisions  are  used.  This  distortion  (Fig.  8A) 
is  well  within  the  limits  which  were  established  for  good  performance 
of  the  PCRK59  element  (see  Section  2) . The  entire  mesh  is  graded 
to  produce  quadrilaterals  with  shapes  as  close  to  square  as  possible 
in  the  bushing  region  (Fig.  8B)  where  the  largest  gradients  in  stress 
and  displacement  are  expected.  Since  the  computer  code  has  been  pro- 
grammed to  confine  the  crack  to  the  bushing,  the  outer  diameter  of 
this  subregion  is  used  to  control  the  potential  extent  of  crack  sizes 
in  analyses  of  single-material  lugs.  The  possible  crack  size  range 
is  given  by: 

a/R1  <_  0.22  for  R2/Rl  1 2 
a/R!  <0.45  for  2.5  < R2//Rl  1 3 
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where 


a = Crack  size 

= Radius  of  bearing  hole 
R^  = Outer  radius  of  lug 

The  stress  intensity  factor  solutions  for  tension  bearing  character- 
istically exhibit  maxima  for  cracks  near  +90°  from  the  lug  axis. 
Figure  8C  illustrates  two  sets  of  solutions:  24  divisions  (coarse 

mesh)  and  32  divisions  (fine  mesh)  around  the  hole.  The  K ^ data, 
plotted  for  0°  £ 0 £ 180°, are  observed  to  have  converged  with  the 
coarse  mesh.  The  KII  data,  plotted  for  180°  £ 0 £ 360°, are  observed 
to  be  near  convergence  with  the  fine  mesh.  Full  plots  are  not  required, 
since  the  solutions  for  tension  bearing  are  always  symmetric  about  the 
lug  axis.  The  upper  part  of  Fig.  8C  also  illustrates  contour  plots 
of  the  stress  components  in  an  uncracked  lug  with  the  same  dimensions. 
These  stress  analyses  were  used  to  demonstrate  a correlation  between 
K^,  for  small  cracks  and  respectively  [8],  and  the 

stress  analysis  option  has  been  retained  in  the  final  program. 

Figure  9 illustrates  the  mesh  details  and  typical  solutions  for 
positive  shear  bearing  (load  position  2 in  Fig.  7)  for  the  same  lug 
dimensions  as  shown  in  Fig.  8B.  Stress  intensity  factor  solutions 
for  both  uniform  and  cosine  bearing  pressures  are  illustrated.  In 
these  cases,  maxima  occur  for  cracks  positioned  at  unexpected 
angles.  Similar  plots  for  several  other  interesting  cases  appear 
in  the  original  report  [8]. 

The  case  of  tension  bearing  was  considered  to  be  of  primary 
interest  in  the  present  work.  Therefore,  an  extended  parameter 
study  was  conducted  for  this  case  to  evaluate  the  sensitivity  of 
stress  intensity  factors  to  a/R^  and  f°r  both  uniform  and 

cosine  bearing  pressure  distributions.  The  numerical  data  obtained 
from  the  lug  program  analyses  were  reduced  to  a form  convenient  for 
design  charts.  Values  of  the  sensitivity  functions  F^(a/R^,  8)  and 
F^I(a/R^,  0)  are  summarized  in  Tables  2 through  5 (uniform  bearing) 
and  6 through  9 (cosine  bearing).  For  specific  applications,  stress 
intensity  factors  may  be  computed  as: 
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KI  - 7ii  Va/V  9>  KII  = Tna  FII<a/Rl'  91 

where  P is  the  total  bearing  force  value.  These  data  were  obtained 
from  lug  models  with  32  divisions  around  the  bearing  hole,  allowing 
the  crack  to  be  positioned  at  11.25°  intervals. 

Figures  10  and  11  illustrate  the  sensitivities  of  for  cracks 
at  0°  and  90°  respectively.  For  cracks  at  0°,  a cosine  pressure 
distribution  is  seen  to  cause  stress  intensities  more  severe  than 
those  caused  by  a uniform  distribution.  The  situation  is  reversed 
for  cracks  at  90°,  where  the  stress  intensities  are  generally  the 
most  severe.  In  both  cases,  the  dimensionless  parameter  R^/R^  for 
the  lug  is  observed  to  have  a significant  effect.  The  most  rapid 
change  in  stress  intensities  is  seen  to  occur  between  1.5  <_  R2//Rl  — 
This  rapid  change  is  attributed  to  finite-width  effects  similar  to 
those  which  have  been  derived  analytically  for  cracks  in  finite-width 
plates  loaded  in  uniform  tension.  The  behavior  of  these  solutions 
should  be  noted  especially  for  R^/R^  = 1.5  and  R^/R^  = 2,  where 
anomalies  occur.  The  data  at  a/R^  = 0.22  are  obtained  from  finite 
element  models  in  which  the  crack-tip  location  begins  to  approach 
close  to  the  opposite  boundary  of  the  PCRK59  element,  a situation 
in  which  significant  errors  are  to  be  expected  (see  Section  2) . 
Therefore,  the  solutions  # for  these  R^/R-^  values  are  considered  to 
be  valid  only  for  a/R^  _<  0.2.  A similar  but  much  less  severe  effect 
can  be  observed  between  0.4  <_  a/R^  £ 0.45  for  R^/R1  = 2.5  and 
R2/R1  = 3. 

The  sensitivity  data  presented  in  Tables  2 through  9 are  also 
plotted  in  Figs.  12  through  27.  It  should  be  noted  that  F^fa/R^,  0) 
assumes  negative  values  at  some  points.  A negative  value  in  the 
linear  elastic  solution  scheme  means  that  the  two  PCRK59  element 
nodes  which  lie  on  opposite  sides  of  the  c rack  surface  have  overlapped 
one  another.  The  proper  interpretation  is  that  a crack  in  such  a 
location  would  be  closed  by  the  applied  loading;  the  actual  negative 
Fj  value  otherwise  has  no  specific  meaning.  Values  for  and  F^ 

may  be  positive  or  negative  according  to  the  sense  of  the  shear  stress 
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field  in  the  vicinity  of  the  crack  tip.  In  this  case,  the  attachment 
lug  program  routinely  computes  absolute  values,  the  sign  being  of  no 
significance  for  interaction  formulas.  Comparison  of  Figs.  12  through 
15  with  Figs.  20  through  23  shows  that  the  closure  effect  for  cracks 
at  0°  appears  only  in  the  case  of  uniform  bearing  for  F^/R^  — 2- 
Cracks  at  180°  are  subject  to  closure  for  both  uniform  and  cosine 
bearing.  However,  180°  cracks  may  be  opened  by  reversed  (compression) 
loading  [8].  The  Kj  maxima  are  seen  to  occur  for  cracks  near  90° 
in  uniform  bearing,  and  for  cracks  near  100°  in  cosine  bearing. 
Comparison  of  K with  KI3.  (e.g.,  Figs.  12  and  16)  indicates  that 
most  of  the  crack  locations  result  in  significant  combined  stress 
intensities,  while  K is  negligible  only  in  very  narrowly  defined 
regions  about  0 = 0°,  90°  and  180°.  When  considering  the  effects 
of  fatigue  loading,  it  has  been  shown  that  a crack  tends  to  propa- 
gate along  a curved  path  when  the  value  is  significant,  the 

curvature  continuing  until  the  crack  is  oriented  perpendicular  to 
a principal  tensile  stress  [17] . Finally,  it  may  be  observed  that 
the  Kj  charts  exhibit  crossovers  at  0=0°  and  180°  when  the  stress 
intensities  are  small  (Figs.  12,  13,  22  and  23).  The  charts 

also  exhibit  crossovers  in  the  range  0°  £ 0 £ 90°  for  small  stress 
intensities.  These  effects  are  attributed  to  the  usual  "noisiness" 
which  is  associated  with  low-stress  regions  in  a finite  element 
solution . 
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Section  IV 


FASTENER  ROW  ANALYSIS 


The  fastener  row  programs  automatically  generate  and  solve 
finite  element  models  of  the  structural  details  shown  in  Fig.  2C 
and  2D.  A typical  mesh  plan  and  example  results  from  the  single 
fastener  row  program  are  illustrated  in  Fig.  28.  The  mesh  plans 
and  solutions  for  double  fastener  rows  are  generally  similar.  The 
single  and  double  fastener  row  programs  make  extensive  use  of  internal 
automatic  substructuring  and  rotation  transformations  to  avoid  repeated 
regeneration  of  similar  elements  during  parametric  analyses.  These 
programs  are  capable  of  analyzing  only  isotropic  material.  The  bear- 
ing loads  in  each  fastener  hole  are  assumed  to  be  cosine-distributed. 
Values  for  the  total  bearing  force  applied  to  each  hole  may  be  chosen 
separately.  Up  to  10  fastener  holes  may  be  accommodated  in  the  single 
row  program,  while  the  double  row  program  permits  up  to  8 holes  per 
row.  One  or  two  cracks  may  be  placed  at  the  damaged  hole,  which  may 
be  chosen  from  any  of  the  fastener  holes  represented  in  the  finite 
element  model.  The  crack  sizes,  hole  radius  and  centerline  spacing 
are  all  input  parameters.  Detailed  documentation  of  the  fastener 
row  programs  and  the  results  of  some  preliminary  accuracy  tests 
were  reported  previously  [10] . The  earlier  demonstration  runs  and 
the  new  results  presented  in  the  present  report  have  all  been  carried 
out  with  models  containing  three  fastener  holes  per  row. 

The  previous  attempt  to  verify  the  accuracy  of  the  fastener 
row  programs  is  considered  to  be  incomplete  because  of  the  lack  of 
any  independent  analytical  solution  with  which  a direct  comparison 
can  be  made.  Therefore,  some  further  indirect  comparisons  were 
conducted  by  studying  the  sensitivity  of  K to  a/R  for  cracks  per- 
pendicular to  the  applied  load  (0=0°  according  to  the  angular  con- 
vention in  these  finite  element  models) . Analytical  solutions  for 
a cracked  fastener  hole  in  an  infinite  medium  are  available  in  the 
forms  [ 3 ] : 

K = o /rFa  F (s) 

I 00 
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for  uniform  remote  tension  loading  o^,  where  s = a/ (R+a) , and: 


K][  = p /rra  F (s) 


for  uniform  pressure  loading  p applied  to  the  periphery  of  the  hole. 
In  order  to  compare  the  finite  element  results  with  these  solutions, 
an  "equivalent  pressure"  corresponding  to  the  distributed  bearing 
load  was  calculated  from  a solution  of  the  biharmonic  equation  for 
an  uncracked  structure,  which  gives: 


°ee (r,e)  = 


2P 

2 

7T  R 


r /Rv  2 2 ,Rx  4 oqi  r 4 /R\  6 2 /Rv  Q i 

[ (-)  + j(-)  COS20  + {jj(-)  - j-g(j)  } cos46  + ...] 


where  R represents  the  fastener  hole  radius  and  r ^ R.  The  "equivalent 
pressure"  was  defined  arbitrarily  for  a crack  at  0 = 0°  by  requiring 
aQQ(r,0)  to  be  equal  to  the  hoop  stress  created  by  a uniform  internal 

U U 

pressure : 


op  ft  P 

aQ0(r,O)  = — 2~  [1+J5+-.-]  - .377  --  = p (equiv.) 

7T  R 

Sensitivity  factors  F(s)  for  the  finite  element  solutions  were  then 
computed  from  the  numerical  K^.  values  by  solution  of  the  equivalent 
relationship : 


Kx  - .377  f-  /Wa  F (s) 

The  comparison,  illustrated  in  Fig.  29,  indicates  that  the  fastener 
hole  loaded  by  a cosine  pressure  distribution  has  a behavior  inter- 
mediate between  the  two  analytical  cases,  but  closer  to  the  case  of 
uniform  internal  pressure  than  to  remote  tension  loading. 

Additional  indirect  verification  was  obtained  by  comparing  the 
finite  element  solution  with  the  case  of  a crack  of  length  of  2a 
in  an  infinite  medium  with  point  loads  + P applied  to  opposite 
surfaces  at  the  center  of  the  crack.  The  analytical  solution  for 
this  case  is  given  by  [3] : 


Kx  = P/Zira 

If  one  plots  a sensitivity  function  F(s)  for  the  case  of  a crack 
emanating  from  a fastener  hole  in  the  form: 

KI  - 7?J  FI<S-  e» 
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then  as  a/R  becomes  large,  it  is  reasonable  to  look  for  asymptotic 
behavior  which  makes  the  stress  intensity  trend  toward: 

KI  * TFTSTa)'  as  a/R  " “ 

The  above  trend  is  equivalent  to: 

Vs'0'  ^ 

Figure  30  illustrates  the  comparison  for  sensitivity  functions 
F^(s,  0)  computed  from  the  values  obtained  in  a series  of  analyses 
with  the  single  fastener  row  program.  The  curve  for  0 = 0°  is  seen 
to  exhibit  the  correct  trend,  but  it  lies  below  the  curve  for  /s. 

The  difference  may  be  attributed  to  the  fact  that  the  cosine-distributed 
applied  loading  is  diffuse  compared  to  the  point  load  considered  in 
the  analytical  solution.  These  comparisons,  taken  together  with  the 
performance  tests  of  the  rectangular  panel  program  (Section  2)  are 
considered  to  be  sufficient  to  justify  confidence  in  the  fastener 
row  program  solutions.  The  rectangular  panel  accuracy  tests  are 
applicable  in  the  present  case  because  the  finitie  element  models 
for  the  fastener  row  details  are  constructed  from  the  same  basic 
building  blocks,  in  the  same  shape  ranges,  as  the  rectangular  panel 
models . 

Another  important  aspect  of  Fig.  30  is  the  crossover  behavior 
which  is  apparent  at  small  crack  sizes  (0.0158  £ s £ .0875).  This 
type  of  behavior  is  even  more  apparent  in  the  K sensitivity  func- 
tions (Fig.  31),  which  also  show  anomalously  high  values  for  F (s,0) 
at  s = .0158.  The  crossovers  in  K cannot  be  considered  to  result 
only  from  normal  solution  "noise"  because  cracks  within  0 = + 30° 
are  involved,  and  these  are  high-stress  regions.  Furthermore,  one 
should  expect  K ^ (and  hence  Fjj)  to  approach  zero  as  the  crack 
size  approaches  zero,  and  this  expectation  is  followed  only  for 
0 = 0°,  -15°  in  Fig.  31.  Hence,  the  solutions  at  s = .0158  are 
considered  to  be  seriously  in  error.  The  cause  is  attributed  to 
the  proximity  of  the  crack-tip  location  to  the  fastener  hole  sur- 
face, leading  to  the  conclusion  previously  mentioned  (Section  2) 
that  the  crack-tip  location  test  of  an  edge-cracked  semi-infinite 
strip  is  an  incomplete  test. 
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Based  on  the  above  considerations,  data  obtained  from  the 
fastener  row  programs  were  considered  to  be  valid  in  the  ranges 
0.05  £ s <_  0.528  for  and  0.1  £ s £ 0.528  for  K^,  and  sensitivity 
charts  were  prepared  accordingly  from  the  results  of  an  extended 
parameter  study.  The  data  were  taken  from  models  with  a 5/8-inch 
fastener  hole  diameter  for  the  seven  crack  sizes  and  corresponding 
s-values  given  in  Table  10.  A 5/8-inch  fastener  diameter  was  selected 
because  this  size  is  found  in  many  different  aircraft.  For  this  fas- 
tener size,  Table  10  indicates  that  the  invalid  K-solutions  (s  = 0.0158) 
correspond  to  crack  sizes  considered  in  the  durability  portion  of  the 
Air  Force  damage  tolerance  criteria,  while  the  lower  limit  for  the 
valid  K-solutions  corresponds  to  crack  sizes  considered  for  struc- 
tural integrity  ( safety-of-f light ) in  the  criteria  [1].  Also,  the 
upper  limit,  which  represents  the  capability  of  the  finite  element 
meshes,  is  seen  to  be  sufficiently  large  to  have  some  confidence 
that  most  of  the  crack  size  range  over  which  the  fastener  hole  has 
a significant  influence  has  been  covered. 

Reduced  data  for  the  sensitivity  functions  F^(s,  0)  and  0) 

for  several  of  the  cases  considered  in  the  extended  parameter  study 
are  given  in  Tables  11  through  17.  In  most  of  the  cases  considered, 
the  fastener  spacing  was  fixed  at  C/D  = 4 (C/R  = 8) . However,  a 
limited  study  of  the  sensitivity  to  C/D  was  conducted  with  the 
single  fastener  row  program  at  values  C/D  = 3.2,  3.6  and  4.8  (Table  12). 
The  applied  loads  were  assigned  equal  values  for  all  cases  except 
those  in  which  the  effect  of  an  unloaded  damaged  hole  was  studied. 

The  data  have  been  reduced  in  a manner  such  that  stress  intensities 
for  specific  cases  can  be  computed  from: 


Ki  = Til  Vs'  0) 


K„ 


FTT(s,  0) 


'II  " TTa  II 

where  P represents  the  bearing  force  applied  to  any  one  fastener 
hole  in  the  structure,  whether  the  damaged  hole  is  loaded  or  unloaded. 

Sensitivity  charts  have  also  been  prepared  from  the  cases 
summarized  in  Tables  11  through  17,  as  well  as  for  several  other 
cases  of  minor  interest  which  have  not  been  tabulated.  The  sensi- 
tivity charts  are  presented  in  Figs.  32  through  58.  For  convenience 
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of  visual  interpolation  in  those  cases  where  solutions  were  obtained 
for  all  seven  crack  sizes,  the  tabulated  data  have  been  used  to  com- 
pute the  data  for  plotting  at  convenient  s-values  (i.e.,  s = 0.05, 

0.1  etc.)  by  linear  interpolation.  All  of  the  solutions  are  obtained 
from  meshes  with  24  divisions  around  the  damaged  hole  (crack  locations 
at  15°  intervals) . 

Figures  32  and  33  illustrate  the  sensitivity  charts  for  a single 
row  of  fasteners  with  the  center  hole  damaged.  In  view  of  the  symmetry 
about  the  vertical  centerline  of  the  structure  and  loading,  these 
charts  cover  only  -90°  <_  0 <_  + 90°.  The  maxima  are  observed  to 
shift  from  0 = -30°  to  0 = -15°  as  the  crack  size  increases.  As 
was  found  in  the  attachment  lug  analyses,  the  values  are  negli- 

gible only  in  fairly  narrow  regions  near  0=0°,+  90°.  Figures 
34  and  35  illustrate  the  sensitivity  for  three  selected  crack  sizes 
as  the  fastener  spacing  C/D  is  varied.  These  charts  have  been  pre- 
pared to  be  used  in  conjunction  with  Figs.  32  and  33.  For  example, 
in  a specific  case  may  be  computed  from: 

KI  = 7ua  FI(s'  6)  Gi(c/D'  s'  0) 

where  F^  is  obtained  from  Fig.  32  and  G ^ is  obtained  from  Fig.  34. 

The  case  C/D  = 4 is  not  plotted  in  Figs.  34  and  35,  since  by 
definition: 

GI(4,  s,  0)  = GII(4,  s,  0)  = 1 

Some  of  the  angular  variations  of  G and  G ^ can  probably  be  attributed 
to  solution  noise.  The  conclusion  of  primary  interest  about  sensitivity 
to  fastener  spacing  concerns  the  behavior  of  for  cracks  in  the  range 
-30°  <_  0 <_  0°.  The  data  in  Fig.  34  are  fairly  consistent  over  this 
range,  indicating  that  the  K value  (in  terms  of  C/D  = 4)  increases 
by  15  to  25  percent  for  C/D  =3.2,  by  7 to  10  percent  for  C/D  = 3.6, 
and  decreases  by  10  to  20  percent  for  C/D  = 4.8. 

Although  a crack  at  a fastener  hole  is  viewed  primarily  from 
the  standpoint  of  its  effect  in  degrading  structural  integrity, 
there  is  a possibility  that  the  presence  of  the  crack  may  introduce 
enough  additional  local  compliance  in  the  structure  to  unload  the 
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fastener.  Therefore,  stress  intensity  factor  solutions  for  rows 
of  fastener  holes  in  which  the  damaged  hole  is  unloaded  are  of  some 
interest.  Figures  36  and  37  illustrate  the  effect  on  and  K^ 
in  a single  row  of  three  fastener  holes  for  three  selected  crack 
sizes.  The  major  point  of  interest  is  that  the  bearing  loads  in 
the  adjacent  holes  appear  to  provide  about  1/6  of  the  total  stress 
intensity  in  the  range  -30°  £ 0 £ 0°.  Thus,  even  a partial 

unloading  of  the  fastener  in  the  damaged  hole  would  have  a dispro- 
portionately beneficial  effect.  One  might  also  conclude  tentatively 
from  these  results  that  the  contributions  of  K^.  from  more  remote 
fastener  holes  tend  to  die  out  rapidly  as  a function  of  the  distance 
between  the  contributor  and  the  damaged  hole.  Before  acceptance, 
this  conclusion  should  be  verified  with  analyses  of  models  with 
five  or  more  holes  per  row. 

Figures  38,  39  and  40  repeat  the  study  of  basic  sensitivities 
for  the  three-fastener  single  row  with  the  left  hole  damaged.  In 
these  cases,  some  asymmetry  appears  in  the  solutions  because  of 
the  location  of  the  damaged  hole  near  a free  lateral  edge.  In 
particular,  the  values  (Fig.  38)  for  cracks  directed  toward  the 
free  edge  (6  near  180°)  are  somewhat  higher  than  the  corresponding 
values  for  interior  cracks  (0  near  0°).  A similar  effect  on  K 
can  be  seen  in  Fig.  39  (compare  0 = 30°  with  150°  and  6 = 330°  with 
210°) . The  sensitivity  of  K to  fastener  spacing  (Fig.  40)  appears 
to  follow  the  behavior  exhibited  in  the  case  of  the  center  hole 
damaged.  (The  data  for  the  run  s = 0.138  were  inadvertantly  lost 
for  0 > 135°.) 

The  remaining  charts  present  the  results  of  a parameter  study 
conducted  with  a double  row  (three  fasteners  per  row)  with  C/D  = 4. 
Sensitivity  to  fastener  spacing  was  not  included  in  this  part  of 
the  study.  Figures  41  and  42  illustrate  a case  in  which  the  center 
hole  in  the  upper  row  is  damaged.  These  solutions  are  similar  to 
the  corresponding  case  of  a single  row  of  fasteners,  except  that 
asymmetry  is  exhibited  because  of  the  staggered  arrangement  of  the 
holes . 
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Figures  43  and  44  illustrate  a case  in  which  the  center  hole 
in  the  lower  row  is  damaged.  Comparison  with  the  previous  case  shows 
that  the  stress  intensity  factors  for  this  case  are  significantly 
higher.  This  happens  because  a crack  in  the  lower  row  is  influenced 
by  the  bearing  loads  in  the  upper  row  as  if  these  loads  were  remote 
tension  [10],  Also,  the  solution  asymmetry  is  much  less  apparent 
in  this  case.  Figures  45  and  46  illustrate  the  effect  on  and  K^ 
when  the  damaged  hole  is  unloaded.  Again,  significant  decreases 
in  stress  intensity  occur,  but  not  to  the  extent  observed  in  the 
case  of  the  single  fastener  row.  Apparently,  the  loads  in  neighbor- 
ing holes  in  the  upper  row  can  contribute  to  the  crack-tip  stress 
intensity  much  more  easily  than  can  the  loads  in  holes  which  are 
separated  laterally  from  the  damage  site. 

The  remaining  charts  illustrate  results  obtained  with  only 
three  selected  crack  sizes  for  some  minor  cases  of  interest.  The 
stress  intensities  for  a damaged  hole  in  the  lower  row  near  the 
right  edge  are  somewhat  lower  than  for  damage  at  the  center  hole, 
and  the  effects  of  unloading  are  similar  (Figs.  47  through  50). 

When  the  left  hole  in  the  upper  row  is  damaged,  the  results  are 
similar  to  the  corresponding  case  of  a single  fastener  row  (Figs. 

38,  39) , but  in  the  present  case  unloading  of  the  damaged  hole 
results  in  large  regions  of  crack  closure  due  to  the  compression 
effects  of  the  bearing  loads  in  the  lower  row  (Figs.  51,  52) . The 
behavior  of  K ^ is  somewhat  anomalous.  Unloading  causes  a large 
decrease  in  the  stress  intensities  for  large  cracks  (s  = 0.528), 
but  a much  smaller  decrease  for  intermediate-sized  cracks  (s  = 0.39). 
as  shown  in  Figs.  53  and  54.  The  final  case,  in  which  the  right  hole 
in  the  upper  row  is  damaged  (Figs.  55  through  58)  is  similar  to  the 
previous  case  if  the  crack  angle  is  reflected  (i.e.,  the  roles  played 
by  -90°  <_  0 <_  + 90°  and  90°  <_  0 £ 270°  are  reversed). 

It  should  be  noted  that  all  of  the  sensitivity  charts  presented 
in  this  report  have  been  plotted  by  making  linear  connections  between 
the  computed  data  points,  unlike  the  fairing  which  was  done  in  the 
previously  reported  polar  plots.  Linear  interpolation  has  been  used 
purposely  in  the  sensitivity  charts  to  avoid  applying  any  bias  to 
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the  intermediate  angular  positions  for  which  the  programs  cannot 
compute  solutions.  Obviously,  the  true  solutions  should  be  smooth 
curves.  In  most  cases,  the  linear  interpolations  are  probably 
close  enough  for  engineering  purposes,  but  it  is  also  evident  that 
the  computed  solutions  have  missed  peaks  in  many  places.  This  is 
particularly  true  for  the  K ^ charts. 

Finally,  it  must  be  remarked  that  the  fastener  row  solutions 
presented  here  probably  do  not  correspond  directly  to  cases  of 
practical  interest,  in  which  the  load  introduction  most  likely 
will  exhibit  some  variation  with  the  location  of  individual  fasteners 
in  the  panel.  However,  the  case  of  equal  loads  was  chosen  for  study 
as  the  best  method  of  illustrating  general  behavior  without  unnec- 
essarily cluttering  the  presentation  with  additional  parameter 
variations.  In  any  specific  case,  the  fastener  row  programs  can 
be  run  with  individually  chosen  loads  at  each  fastener  location. 
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Section  V 


DISCUSSION  AND  CONCLUSIONS 

The  results  of  a research  program  involving  application  of 
the  assumed-stress  hybrid  finite  element  method  to  the  analysis 
of  stress  intensity  factors  in  some  typical  aircraft  structural 
details  have  been  presented.  The  key  feature  of  the  hybrid  method 
which  has  made  this  work  possible  is  a special-purpose  element 
which  contains  the  crack  tip  and  its  stress  singularity,  but  which 
may  be  coupled  to  conventional  finite  elements  in  a standard  Matrix 
Displacement  Method  analysis  scheme.  A secondary  feature,  which 
was  developed  in  the  course  of  the  present  program,  is  another 
special-purpose  hybrid  element  which  permits  the  rapid  grading 
of  a mesh  from  a very  coarse  cartesian  pattern  to  a locally  fine 
polar  pattern  for  surrounding  fastener  holes.  Tailored  modular 
computer  codes  were  prepared  with  these  building  blocks  to  create 
automatically  generated  finite  element  models  of  four  classes  of 
structural  details:  an  attachment  lug  loaded  in  bearing  through 

its  pin  hole;  a rectangular  stiffened  or  unstiffened  panel  loaded 
by  remote  tension  containing  a centered  or  offset  open  fastener 
hole;  a single  row  of  fastener  holes  loaded  in  bearing;  and  a 
double  staggered  row  of  fastener  holes  loaded  in  bearing. 

Extended  parameter  studies  were  carried  out  with  these 
computer  codes  to  assess  the  general  behavior  of  the  stress  in- 
tensity factor  solutions  and  their  sensitivity  to  crack  size, 
fastener  spacing,  etc.  The  results  of  the  parameter  studies  have 
been  presented  in  the  form  of  handbook  charts  which  can  be  used 
to  compute  stress  intensity  factors  for  specific  cases  of  interest. 
Solution  accuracy  was  verified,  where  possible,  by  comparison  with 
independent  analytical  solutions.  In  those  cases  for  which  ana- 
lytical solutions  are  not  yet  available,  verification  was  accom- 
plished by  tracking  the  performance  of  individual  finite  elements 
over  ranges  of  distortions  in  shape  and  size  which  covered  the 
ranges  used  in  the  computer  codes.  Also,  experience  with  the 


22 


codes  during  the  parameter  studies  demonstrated  that  solutions 
could  be  obtained  for  a complete  parametric  variation  of  the 
angular  position  of  a crack  with  moderate  core  storage  require- 
ments and  execution  times  of  the  order  of  one  minute.  The  com- 
puter codes  are  therefore  considered  to  be  economically  practical 
analysis  tools  which  give  reasonable  answers  for  engineering 
purposes. 

Some  limitations  of  the  finite  element  building  blocks  and 
in  the  general  computer  codes  were  encountered  in  the  course  of 
the  extended  parameter  studies.  The  hybrid  crack  element  appears 
to  give  good  solutions  when  the  crack-tip  location  is  restricted 
to  a range  between  25  and  75  percent  of  the  total  distance  across 
the  crack  element.  Other  limitations,  discussed  in  detail  below, 
also  lead  to  some  ideas  for  future  improvements. 

The  attachment  lug  computer  code  was  found  to  provide  well- 
converged  solutions  over  the  range  of  interest  in  the  ratio  of 
lug  width  to  bearing  hole  diameter.  However,  the  present  code 
is  restricted  somewhat  in  the  extent  of  crack  sizes  which  may  be 
analyzed.  This  restriction  arises  only  from  the  way  in  which  the 
mesh  generation  scheme  was  programmed,  not  from  any  fundamental 
difficulty  with  the  analysis  method.  A routine  reprogramming 
of  the  mesh  generation  scheme  would  permit  the  study  of  cracks 
which  approach  the  outer  boundary  of  the  lug. 

The  solutions  obtained  from  the  rectangular  panel  computer 
code  were  found  to  be  uninteresting,  in  that  the  presence  of 
stiffeners  and/or  an  offset  of  the  fastener  hole  had  very  little 
effect  on  the  stress  intensity  factors.  The  solutions  for  a 
centered  hole  in  an  unstiffened  panel  were  found  to  be  quite 
useful  in  assessing  the  performance  of  the  special-purpose  hybrid 
elements.  However,  it  has  been  concluded  that  too  much  is  demanded 
of  a mesh  generation  scheme  which  attempts  to  model  a crack  at  a 
fastener  hole  together  with  the  crack  tip  closely  approaching  a 
stiffener  or  an  unstiffened  free  edge.  These  latter  cases  appear 
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to  be  of  much  greater  interest,  and  future  work  in  this  area 
should  concentrate  on  a small  subregion  of  the  panel,  with  the 
crack  location  restricted  to  little  or  no  angular  variation  about 
the  perpendicular  to  the  applied  loading.  Also,  the  use  of  a 
bi-material  crack  element  should  be  considered  for  crack  tips 
which  lie  close  to  stiffeners.  An  element  of  this  type  has 
already  been  developed  in  another  related  program  of  research. 

The  performance  of  the  special-purpose  fastener  hole  ele- 
ment was  tested  extensively.  It  was  found  to  give  accurate 
answers  down  to  a lower  limit  of  1.27  in  the  ratio  of  its  outer 
edge  dimension  to  its  inner  diameter.  With  four  rings  of  con- 
ventional quadrilaterals  and  the  crack  element  placed  inside  to 
model  the  material  immediately  around  the  fastener  hole,  this 
results  in  a lower  limit  of  3.2  in  the  ratio  of  fastener  spacing 
to  fastener  diameter.  The  special  element  performance  tests  also 
indicated  that  the  spacing  ratio  could  be  pushed  down  to  3.0, 
probably  with  some  loss  in  solution  accuracy.  The  primary  im- 
pact of  these  limits  occurs  in  the  fastener  row  computer  codes, 
which  are  intended  to  analyze  chordwise  wing  skin  splices. 

While  a spacing  ratio  of  3.0  to  3.2  is  probably  adequate 
for  many  cases,  very  highly  stressed  splices  such  as  are  found 
in  the  wing  roots  of  the  Boeing  747  and  Lockheed  C-5A  are  designed 
with  smaller  spacing  ratios  and  cannot  be  analyzed  with  the  com- 
puter codes  in  their  present  forms.  An  improvement  of  the  codes 
in  this  respect  would  appear  to  be  desirable,  since  highly  stressed 
chordwise  splices  are  likely  to  continue  to  be  of  great  interest 
in  future  considerations  of  airframe  structural  integrity.  The 
most  straightforward  way  to  gain  this  improvement  would  be  to 
reduce  the  number  of  quadrilateral  rings  in  the  interior  mesh 
from  four  to  two  (at  least  two  rings  are  required  to  accommodate 
the  crack  element) . If  the  approach  which  was  adopted  for  pro- 
portioning the  interior  mesh  is  not  changed,  the  reduction  from 
four  to  two  rings  would  lead  to  a capability  to  analyze  details 
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with  fastener  spacing  ratios  as  small  as  1.7  to  1.8,  without 
requring  any  modifications  to  the  special-purpose  hole  element. 
However,  a new  series  of  verification  tests  using  the  rectangular 
panel  program  should  be  conducted  if  this  approach  is  adopted. 
Also,  it  should  be  noted  that  a trade-off  in  capabilities  would 
result.  With  the  present  four-ring  interior  mesh,  it  is  possible 
to  analyze  cracks  as  long  as  the  fastener  hole  radius;  if  the 
interior  mesh  were  reduced  to  two  rings,  the  maximum  crack  length 
would  be  reduced  to  about  30  percent  of  the  fastener  hole  radius 
if  no  other  modifications  were  made  in  the  existing  codes. 

Finally,  some  remarks  are  in  order  regarding  the  uses  to 
which  the  solutions  presented  in  this  report  may  be  put.  The 
application  of  fracture  mechanics  to  airframe  damage  tolerance 
analysis  follows  two  main  avenues:  the  prediction  of  failure 

loads  based  on  comparisons  of  stress  intensity  factors  with 
material  fracture  toughness  properties  and  the  prediction  a 
crack  sizes  as  functions  of  time  during  stable  crack  propagation 
under  fatigue  loading.  In  both  cases,  the  stress  intensity 
factor  solutions  of  linear  elastic  fracture  mechanics  are  asso- 
ciated with  an  implicit  assumption  that  the  next  increment  of 
stable  or  unstable  crack  propagation  will  occur  tangent  to  the 
orientation  of  the  crack  at  its  tip.  However,  crack  propagation 
experiments  have  shown  that  fatigue  cracks  tend  to  follow  curved 
paths  in  regions  where  they  are  subjected  simultaneously  to 
mode  I and  mode  II  stress  intensities,  the  curvature  continuing 
until  the  crack  is  advancing  in  a direction  perpendicular  to  a 
principal  tensile  stress.  Therefore,  some  caution  should  be 
exercised  in  applying  the  results  presented  in  this  report  to 
fatigue  crack  growth  analysis.  It  so  happens  that  for  most  of 
the  cases  considered,  the  worst  or  nearly  worst  situation  of 
mode  I stress  intensity  coincides  with  negligible  mode  II  stress 
intensity.  Hence,  cracks  which  are  assumed  to  begin  in  these 
angular  positions  may  also  be  safely  assumed  to  propagate  in 
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straight  lines,  and  crack  growth  calculations  can  be  conducted. 
However,  it  is  not  valid  to  use  the  stress  intensity  solutions 
by  taking  values  for  a series  of  cracks  of  increasing  length 
and  varying  angular  position  in  an  attempt  to  trace  the  tip  of 
a crack  which  propagates  along  a curved  path.  The  restriction 
on  use  of  the  solutions  for  failure  load  prediction  is  somewhat 
less  severe,  since  interaction  formulas  may  be  used  for  straight 
initial  cracks  which  occupy  angular  positions  such  that  both 
mode  I and  mode  II  stress  intensities  are  present. 


26 


Figure  1.  Portion  of  Finite  Element  Model  Surrounding  a Crack  Tip 
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Figure  3 


0.2  0.4  0.6  0.8  1.0 

a/b 


Error  in  Computed  Stress  Intensity  Factor  as  a Function  of 
Crack-Tip  Position 
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PROBLEM 
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(C)  EXAMPLE  SOLUTION 

Figure  4.  Mesh  Construction  and  Example  Solution  for  Rectangular 
Panel 
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Figure  5.  Variation  of  K Solution  Accuracy  with  Crack  Size  (Hole  Element 
Shape  W/D  Fixed) 


Error 


STRESSES  IN  UNCRACKED  LUG 
(1.0  K si  CONTOUR  INTERVALS) 
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Example  Results  for  Shear  Bearing 


(a/R  ,6)  at  0=90 


(a/R  ,0) 


0 (Deg) 


Figure  12.  Attachment  Lug  KJ  Chart  (R2/R;l=1.5,  Uniform  Bearing) 


(a/R  , 9) 


0 (Deg) 


Figure  13.  Attachment  Lug  K Chart  (R  /R  =2,  Uniform  Bearing) 


0 (Deg) 


Figure  14.  Attachment  Lug  Chart  (R  /R  =2.5,  Uniform  Bearing) 
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0 45  90  135  180 


0 (Deg) 

Figure  16.  Attachment  Lug  K Chart  (R  /R  =1.5,  Uniform  Bearing) 


0 (Deg) 

Figure  17.  Attachment  Lug  K Chart  (R  /R  =2,  Uniform  Bearing) 


(a/R 


0 45  90  135  180 

0 (Deg) 


Figure  18.  Attachment  Lug  K Chart  (R  /R  =2.5,  Uniform  Bearing) 


Figure  19.  Attachment  Lug  Chart  (R2/R^=3,  Uniform  Bearing) 


0 (Deg) 


Figure  20.  Attachment  Lug  K Chart  (R  /R  =1.5,  Cosine  Bearing) 

X ^ 1 


0 (Deg) 


Figure  21.  Attachment  Lug  K Chart  (R  /R  =2,  Cosine  Bearing) 


0 (Deg) 

Figure  22.  Attachment  Lug  K Chart  (R/R  =2.5,  Cosine  Bearing) 


0 (Deg) 


Figure  23.  Attachment  Lug  K Chart  (R  /R  =3,  Cosine  Bearing) 
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Figure  24.  Attachment  Lug  K Chart  (R  /R  =1.5,  Cosine  Bearing) 


0 45  90  135  180 

0 (Deg) 


Figure  26.  Attachment  Lug  Chart  (R  /R  =2.5,  Cosine  Bearing) 


0 (Deer) 


Figure  27.  Attachment  Lug  K Chart  (R  /R  =3,  Cosine  Bearing) 


(A)  TYPICAL  FINITE  ELEMENT  MESH 


KSl/Ttf  KSiyitf 

(B)  EXAMPLE  RESULT  (CENTER  HOLE  (C)  EXAMPLE  RESULT  (RIGHT  HOLE 

DAMAGED)  DAMAGED) 


Figure  28.  Single  Fastener  Row  Program  Mesh  and  Example  Results 
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(S)  A 


Figure  29.  Comparison  of  a/R  Sensitivity  with  Analytical  Solutions 
for  Similar  Problems 
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II 


0 (Deg) 


Figure  35.  Sensitivity  of  K. 


(0'S) 


(0's) 


Figure  47.  Double  Row  K];  Chart  (Lower  Right  Hole  Damaged,  C/D=4) 
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.2 


6 (Deg) 


Figure  48.  Double  Row  K Chart  (Lower  Right  Hole  Damaged  and 
Unloaded,  C/D-4) 
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Figure  49.  Double  Row  K Chart  (Lower  Right  Hole  Damaged,  C/D=4) 


0 (Deg) 


Figure  50.  Double  Row  K Chart  (Lower  Right  Hole  Damaged  and 
Unloaded,  C/D±4) 


E 


.8 


0 (Deg) 

Figure  57.  Double  Row  K Chart  (Upper  Right  Hole  Damaged,  C/D=4) 


0 (Deg) 


Figure  58.  Double  Row  K Chart  (Upper  Right  Hole  Damaged  and 
Unloaded,  C/D=4) 


TABLE  1 


TYPICAL  CORE  STORAGE  AND  CPU  TIME  REQUIREMENTS 


Program 

Core  Storage  ^ 

(1/2) 

CPU 

Time  (Min.) 

KBYTES 

Decimal 

Words 

Octal 

Words 

Attachment 

250 

62,500 

172,044 

1.1 

Lug 

Rectangular 

260 

65,000 

176,750 

0.8 

Panel 

to 

to 

to 

to 

300 

75,000 

222,370 

1.0 

Single 

0.25 

Fastener 

272 

68,000 

204,640 

to 

Row  ( 3 ) 

! 

1.25 

Double 

0.25 

Fastener 

344 

86,000 

247,760 

to 

Row  ( 3 ) 

1 

l 

i 

1.25 

(1)  Includes  data  and  object  code  produced  by  IBM 
FORTRAN  G-l  and  H(0)  compilers. 

(2)  Times  are  quoted  for  one  complete  set  of  K^,  K__, 
solutions,  i.e.,  one  set  of  dimensions  and'L 

all  angular  crack  positions,  on  an  IBM  S-370/168. 

(3)  Time  and  storage  quoted  for  models  possessing  three 
fastener  holes  per  row. 
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TABLE  2 


UNIFORM  BEARING  (R2/R  =1.5) 


.02 

.04 

.05 

.10 

.15 

. 20 

.22 

0 

-.069 

0 

-.089 

0 

-.105 

0 

-.130 

0 

-.082 

0 

-.006 

.001 

■EBi 

| 

11.25 

.006 

.011 

.022 

.040 

.022 

.047 

.050 

.109 

.115 

.168 

.168 

.210 

.811 

.314 

22.5 

■ 

.228 

.028 

.255 

.031 

. 377 
.099 

.473 

.184 

.495 

.277 

1.004 

.482 

33.75 

.227 

.059 

.342 

.016 

.384 

.019 

.560 

.047 

.682 

.161 

.706 
. 326 

1.046 

.626 

45 

.241 

.079 

. 358 
. 026 

.404 

.021 

.610 

.090 

.760 

.257 

. 815 
.512 

.979 

.888 

56.25 

.244 

.072 

. 372 
.018 

.428 

.042 

.683 

.231 

.884 

.461 

1.005 
. 799 

.988 

1.212 

67.5 

.256 

.049 

.411 

.058 

. 487 
.095 

.834 
. 320 

1.116 

.571 

1. 305 
.918 

1.129 
1.  303 

78.75 

.234 

.024 

.404 

.043 

.487 

.068 

. 885 
.207 

1.226 

.357 

1.495 

.566 

1.228 

.793 

90 

.196 

.006 

.348 
. Oil 

.428 

.015 

. 801 
. 045 

1.137 

.088 

1.424 

.136 

1.162 

.124 

101.25 

.167 

.006 

. 288 
.064 

. 354 
.095 

.666 

.262 

. 938 
.447 

1.163 
. 679 

.921 
. 793 

112.5 

.121 

.012 

.193 

.077 

.235 

.114 

. 435 
. 309 

.606 

.514 

m 

.560 
. 888 

123.75 

.062 

.004 

.089 

.044 

.109 

.072 

.200 
. 224 

.278 
. 382 

.348 
. 570 

.212 
. 657 

135 

.024 

.010 

.028 

.008 

.034 

.026 

.056 

.134 

.071 
. 249 

.089 
. 380 

-.017 

.434 

146.25 

.001 

.006 

-.009 

.013 

-.012 

.031 

BjftEQ 

mm 

-.064 

.227 

-.079 

.340 

-.156 

.387 

157.5 

-.024 

.002 

m 

-.078 

.039 

-.149 

.123 

-.212 

.202 

-.252 

.292 

-.301 

.330 

168.75 

-.058 

0 

-.106 

.018 

-.129 
. 029 

-.236 

.080 

-.322 

.125 

-.376 

.175 

-.404 

.198 

180 

-.066 

0 

-.121 

0 

-.268 

0 

-.360 

0 

-.420 

0 

-.440 

0 

K = — F(a/R, ,6)  Upper  values  - F 

/ 1 

>/1Ta  Lower  values  - F^ 
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UNIFORM 


= — F (a/R, / e ) 

/ — X 

/it  a 


K 


TABLE  3 


BEARING  (R2/R1=2) 


.05 

.10 

.15 

.20 

.22 

.036 

-.043 

-.030 

-.034 

0 

0 

0 

0 

■EM 

.009 

-.007 

mm 

-.006 

. 312 

.037 

.069 

■SB 

.104 

.139 

.051 

.073 

.086 

.063 

.341 

.044 

.095 

.142 

.179 

.250 

.107 

.152 

.174 

.154 

.375 

.026 

.166 

. 240 

.344 

.141 

.216 

.264 

.270 

.422 

.016 

.104 

.206 

.320 

.441 

.167 

.281 

.371 

.421 

.497 

.027 

.140 

.260 

.399 

. 518 

.194 

.352 

.484 

.577 

.574 

.043 

.156 

.268 

.395 

.493 

.215 

.402 

.558 

.677 

. 609 

.034 

.104 

.168 

.238 

.299 

.230 

.419 

.573 

.695 

. 586 

.006 

.014 

.025 

.038 

.019 

’.YlT 

. 386 

.523 

.627' 

II 

.050 

.130 

.205 

.288 

II 

.179 

.309 

. 412 

.490 

.369 

.069 

.178 

.277 

. 386 

.403 

7123” 

.266 

7 315 

. 207 

.057 

.155 

.247 

.346 

.358 

.069 

.134 

.157 

.065 

.034 

.106 

.177 

.255 

. 258 

7021 

.024 

023 

' .026' 

-.048 

.028 

.086 

.144 

.208 

.212 

.027 

-.055 

-.143 

.026 

.073 

.119 

.169  | 

.172 

.064 

-.115 

-.154 

,-.171 

KnSI 

.017 

.046 

.072 

.100 

.103 

.078 

-.138 

-.182 

-.200 

-.233 

0 

0 

0 

0 

0 

Upper  values  - F^ 
Lower  values  - F^ 
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TABLE  4 


UNIFORM  BEARING  (R2/R  =2.5) 


\a/R1 
0° \ 

CM 

o 

. 

o 

o 

Ln 

.10 

.15 

. 20 

. 25 

. 30 

. 35 

.40 

| .45 

0 

-.015 

0 

-.017 

0 

-.019 

0 

-.015 

0 

.001 

0 

.022 

0 

.047 

0 

.075 

0 

.104 

0 

.129 

0 

.182 

.001 

11.25 

-.013 

.010 

-.009 

.038 

-.009 

.043 

.004 

.074 

.025 

.099 

.047 

.113 

Tibi 

.138 

.069 

.118 

.092 

.113 

.118 

.103 

141 

.088 

.200 

.076 

22.5 

-.004 

.002 

.012 

.032 

.002 

.034 

.048 

.069 

.077 

.108 
_j 

.119 

.158 

.133 

.167 

. 150 
.169 

.171 

.166 

.244 

.175 

33.75 

.017 

.044 

.052 

. 095 

.134 

. 162 

.181 

.193 

. 208 

. 230 

. 309 

.037 

.015 

.018 

.011 

.058 

.102 

.141 

.176 

.211 

. 249 

.309 

45 

.044 

.08  0^ 

.091 

. 147 

. 199 1 

. 240 

. 270  1 

.293 

. 316 

.343 

.415 

.076 

.063 

.065 

.021 

. 038 

.090 

.142 

.198 

. 266 

. 347 

.465 

56  25 

.064 

.107 

.123 

.206 

.282 

. 346 

. 400 

.444 

.483 

.517 

.569 

.095 

.079 

.073 

.002 

.074 

.129 

.183 

.248 

. 331 

.436 

.580 

67.5 

.067 

.117 

.139 

.253 

. 356 

. 444 

. 522 

.590 

.649 

. 696 

. 72  7 

.081 

,057 

.045 

.043 

. 110 

.157 

.201 

..  ...J 

. 256 
' .665’ 

.330 

.422 

. 548 

78.75 

.059 

^ .115 

.14  3' 



. 281  I 

. 397 

.496 

.585 

. 739  1 

.797* 

. 817 

. 042 

,022 

- i 

.011 

. 049' 

j 

.088 

.111 

.132 

...  ■ 

.159 

.196 

. 244 

. 310 

90 

.067 

. 128 



.160 

. 30  7, 

i 4 23 

' .518  i 

. 604  1 

" . 6 84  ! 

. 758  t 

” .*819” 

~ .833" 

.002 

.002 

. 003 

. 007: 

.013 

.019  I 

. 027  j 

. 034  1 

j 

.042  : 

.050 

.056 

; .0  8 5 

; . 14  3; 

.174; 

i .3161 

.4  28 

.520  j 

. 599  | 

. 668  I 

' . 728  ' 

” . 7 73  ' 

. 7 72  “ 

101. 25 

.022 

.001! 

1 1 

. 0 1 6 j 

i -088j 

.136 

. 169  ! 

.......  i 

. 200  | 

. 2 35  j 

. 278  ! 

.329  j . 390  j 

112.5 

1 “.’lib 

.160. 

'.188': 

.304 

i 39  8 

.475  j 

. 5 39  f 

’ .591' 

"Ten”; 

. 655 

.642 

.026 

. 002  [ 

. 0 24  ! 

. 129  ’ 

i 

.201 

.250 

. 293  ] 

. 342  1 

.401  | 

. 474  , 

. 560  j 

.ii2 

. 1 50  | 

:i7i! 

.2511 

.315;  .3681 

.410 ; 

' . 44  2 ; 

. 46  2 j 

. 469  ; 

.448  ! 

12  3. 7 5 

. 01 8 : 

,001  j 

. 019  i 

- - | 

. 1 1 3 j 

.184  j 

.234 

. 279  j 

. 329  j 

. 391 

. 463  | 

.545 

135 

.084 

.108 

.1201 

.167] 

. 201  ! 

.230 

. 253 

. 266  I 

.| 

. 271  ! 

. 269  ' 

.247 

.008 

,002  ! 

i 

.011' 

i 

. 073] 

.125  ! 

.168  j 

.208] 

. 253  | 

..  4 

. 307  j 

. 366  - 

.432 

” . 041" 

. 0 50  j 

. 055 

. 071 

. 083  ! 

. 092  ; 

. 098  ! 

. 099  i 

. 095  | 

. 089  ; 

. 074  ; 

14  6.25i 

.005 

.004  j 

.00  5] 

.0  49 

. 088  ; 

.121  f 

. 155  ! 

.192 

.235  ; 

. 283  j 

. 335  j 

157.5 

.0051 

-.001  ! 

-.  002' 

-.014 

-- 022 i 

-.029 



-.  035  j 

-.042  ' 

-.  050  : 

-.054  ; 

-.060  ! 

.008 

.004  , 
-.030 J 

.004! 

.041! 

.071 ! 

. 096  ; 

.121  | 

.148  j 

.179  ! 

.212 

. 248  j 

168 . 75 j 

-".015 

- . 0 3 8 ! 

-.0  71* 

-.095  J 

-.1121 

-.127  | 

-.139  i 

-.148  ! 

-.152 

-.151  | 

.007 

, 003 ; 

.003; 

. 027; 

. 045  j 

.059 

.072  j 
- . 159  1 

. 086  ] 

.102  j 

. 119 

.137  [ 

1 

- ! 0 21 

-.040  I 

-.049' 

- • 09 li 

-•121 j 

-.142” 

-.173  j 

183  i 

-.186] 

-.183”] 

180 

0 

__°J 

0 

0 I 

o ! 

. i 

0 

• i 

» J 

. 0 J 

. 0 J 

...  °j 

K = — F(a/R  ,0)  Upper  values  - F_ 

1 T . * 

Lower  values  - F 
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TABLE  6 


COSINE  BEARING  (R  /R  =1.5) 


6 \ 

.02 

.04 

.05 

.10 

.15 

.20 

. 22 

0 

. 037 
0 

.092 

0 

. 116 
0 

. 268 
0 

.458 

0 

. 639 
0 

1.038 

0 

11.25 

. 079 
.004 

.150 
. 008 

.180 

.001 

. 347 
.022 

. 528 
.058 

.677 

.137 

1.062 

.139 

22.5 

.156 

.011 

. 255 
. 003 

. 293 
.011 

.479 

.051 

.639 

.094 

. 731 
.184 

1.079 

.148 

33.75 

.199 

.026 

.309 

.014 

. 349 
.027 

. 534 
.053 

. 675 
.066 

.736 

.095 

1.038 

.005 

45 

.197 

.040 

. 303 
.018 

.341 

.023 

. 519 
. 003 

.654 

.038 

.714 

.102 

.9  46 
. 265 

56.25 

.172 
. 042 

. 273 
.004 

.312 

.004 

. 500 
.080 

.651 

.178 

. 742 
.333 

. 888 
.545 

67.5 

.137 

.030 

.236 

.013 

.279 

.029 

.497 

.133 

. 692 
.264 

.846 

.466 

.9  05 
.686 

78.75 

.101 

.021 

. 200 
.007 

. 249 
.020 

. 505 
.101 

.754 

.207 

.981 

.365 

.971 
. 518 

90 

.131 

.008 

. 243 
.008 

. 301 
.008 

. 588 
.006 

. 863 
.001 

f.115 

.001 

1.046 

.025 

101.25 

.214 

.012 

""  .'346” 
.057 

.4 16 
.082 

" . 722 
.199 

~ . 9 86 
. 317 

1.179 

.472 

1.062 

.539 

112.5 

.2  34 
.017 

.353' 

.082 

.416 

.120 

.6  83 
.299 

78  84” 
.470 

1.005 

.670 

. 880 
.751 

123.75 

. 182 
.025 

. 258 
. 074 

.299 

.106 

.470 

.254 

.597 

.390 

.669 

.546 

. 557 
.598 

135 

.116 
. 016  j 

.157 

.038 

.179 

.057 

. 271 
.146 

.337 

.235 

.373 
. 336 

. 279 
. 361 

146.25 

.0621 

.009 

.080 
. 020 

’ .09  O’ 
.031 

.129 

.090 

.154 

.151 

. 223 

.093 

.239 

157.5  | 

.015 

.005 

.012 

.015 

.012 
. 024 

.008 

.069 

.003 

.116 

.004 

.170 

-.051 

.183 

168.75 

| 

-".019 

.001 

-.038 

.009 

-.04  6 
.015 

-.0  82 
.042 

-.108" 

.069 

-TfiT 

.099 

-.154 

.108 

180 

-.031 

0 

- .055 
0 

-.066 

0 

- . 114 
0 

-.14  7 
0 

-.157 

0 

-.190 

0 

p 

K = — F(a/R..,0)  Upper  values  - F 

-L 


Lower  values  - F 
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TABLE  7 


COSINE  BEARING  (R2/R  =2) 


.02 

.04 

.05 

.10 

.15 

.20 

.22 

0 

.022 

0 

.055 

.068 

0 

.143 

0 

. 219 
_ 0 

.263 

0 

.521 

0 

11.25 

.030 

.002 

.067 

.014 

.081 

.011 

.156 

.003 

.226 

.008 

.263 

.033 

.516 

.019 

22.5 

.048 

0 

.093 

.017 

.108 

.010 

.181 

.002 

.238 

.010 

.256 

.038 

.496 

.005 

33.75 

.061 

.009 

.111 

.009 

.125 

.004 

.194 

.005 

.238 

.009 

. 

.241 

.004 

.460 

.060 

45 

.061 

.020 

.110 

.001 

.124 

.001 

.190 

.021 

.234 

.050 

.237 

.082 

.427 

.156 

56.25 

„ - 

.049 

.025 

.096 

.005 

.111 

.001 

.187 

.042 

.244 

.093 

.270 

.157 

.421 

.238 

67.5 

. . 

.037 

.019 

.083 

.003 

.101 

.005 

.196 

.052 

.281 

.108 

.343 

.180 

.450 

.256 

78.75 

.034 

.017 

.085 

.007 

.109 

.002 

.234 

.035 

.348 

.082 

.441 

.138 

.491 

.195 

90 

.070 

.008 

.132 

.008 

.163 

.008 

.310 

.004 

.435 

.006 

.536 

.015 

.529 
! .033 

101.25 

"Tier 

.010 

.186 

.034 

.222 

.047 

376 

.102 

”7498 

.146 

’.’586’ 

.195 

j~ . 5 3 6~” 

.200 

112.5 

.130 

.008 

| 

.194 

.041 

'.221 

.062 

.156 

.480 

.235 

.548' 

.318 

".475 

.335 

123.75 

.118 

.015 

’ ".‘165 
.040 

.190 

.059 

’.298’ 

.143 

7379 

.216 

______ 

.294 

”7  3 49“ 

.305 

135 

089 

.014 

.119 

.024 

’ .135 
.037 

To  6 

.091 

7258’ 

.142 

’ 77  9 o' 

.199 

.210 

.200 

146.25 

" .054 
.011 

.070 

.014 

.079  ’ 
.021 

'.’117' 

.055 

" 7l45 

.090 

7l65 

.131 

.090 

.129 

157.5 

.020 

.007 

i ~ .023 
.009 

.026 

.014 

.037 

.038 

7 6 46 

.064 

7 058 
.093 

-.008 

.092 

168.75 

-7005 

.002 

-~.012 

.004 

-.014 

.008 

-.024 

.022 

-76 2 8 

.037 

-7  019 
.053 

-.076 

.053 

180 

-.013 

0 

-.025 

0 

-.029 



-.046 

0 

-.055 

0 

-.047 

0 

-.100 

0 

K = — F (a/R,  ,6) 

/FI  1 
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Upper  values  - Fj 
Lower  values  - F 


0 

11.25 

22.5 

33.75 

45 

56.25 

67.5 

78.75 

! 

90 

101.25 

112.5 

123.75 
135 

146.25 

157.5 

168.75 
30 


TABLE  9 


COSINE  BEARING  (R2/R  =3) 


.02 

.04 

.05 

.10 

.15 

. 20 

.25 

. 30 

.35 

.40 

. 45 

0 

.014 

0 

.036 

0 

.048 

0 

.107 

0 

.209 

0 

.253 

0 

.295 

0 

.333 

0 

11.25 

. 015 
. 021 

.039 
. 035 

.050 

.033 

.109 
. 030 

.163 
. 031 

.209 

.030 

.250 

.023 

.285 

.007 

.319 

.015 

.344 

.043 

■ 

22.5 

. 015 
.023 

.042 

.041 

.052 

.036 

.109 

.025 

.160 

.025 

.199 

.027 

.230 

.021 

. 254 
.005 

.279 

.022 

.302 

.057 

Wm 

33.75 

.013 

.003 

.040 

.016 

.049 

.009 

.096 

.006 

.139 

.004 

.171 

.003 

. 193 
.007 

.209 

.003 

. 226 
.008 

.321 

.035 

45 

.007 

.028 

.032 

.023 

.038 

.029 

.078 

.032 

.116 

.023 

.146 

.013 

.166 

.002 

.181 

.010 

.197 

.025 



.224 

.040 

.305 

.068 

56.25 

.021 
. 053 

.027 
* 056 

.070 
. 045 

.114 
. 032 

.150 
. 020 

.178 
. 0Q4 

.201 
. 022 

. 227 
. 058 

.259 
• 104 

.335 
. 164 

67.5 

-.005 

.052 

.016 

.057 

.024 

.058 

.135 

.032 

.184 

.022 

.229 

.002 

.268 
. 031 

.309 

.078 

. 348 
.136 

78,75 

-.005 

.046 

.019 

.052 

.031 

.052 

.107 

.040 

.181 

.027 

. 246 
.012 

.306 

.011 

. 359 
.043 

.410 

.082 

.449 

.128 

.479 

.181 

90 

.035 

.009 

.071 

.011 

.090 

.012 

.184 

.010 

. 262 
.004 

.328 
. 006 

.387 
. 017 

.442 

.028 

.490 

.038 

.528 

.044 

.539 

.051 

101.25 

.087 

.001 

.136 

.017 

’’  161 
.030 

.267 

.078 

.345 

.100 

. 407 
.107 

. 459 
.111 

. 504 
.116 

. 54  2 
.127 

. 568 
.144 

. 565 
.165 

112.5 

.114 

.017 

.164 

_jl006_ 

.189 

.025 

. 295 
.108 

.373 

.159 

.435 

.187 

.483 

.208 

.519 

.230 

.542 

^258 

.550 

.293 

.527 

.331 

123.75 

.127 

.005 

.171 

.013 

.192 

.031 

. 387 
.189 

.425 

.214 

.451 
. 240 

.462 

.272 

.460 

.309 

.431 

.348 

135 

HH 

.141" 

.016 

.157 
. 028 

.216 

.074 

Bl 

'T293- 

.130 

1 ..  . 

.319 

.152 

"73  35 
.177 

~340 
. 205 

.335 

.235 

.308 

.265 

146.25 

.069 

.013 

~T088~ 

.012 

.098 

.018 

.135 

.040 

.163 

.056 

.186 

.072 

KH 

.214 

.111 

.216 

.133 

. 214 
.156 

.194 

.175 

157.5 

m 

.038 

.004 

.042 

.008 

.060 

.019 

.075 

.029 

.089 

.041 

.101 

.056 

.109 

.070 

.113 

.085 

.115 

.100 

.108 

.110 

168.75 

.004 

.001 

.004 

.001 

.005 

.001 

.014 

.014 

.021 

.022 

.028 

.030 

.035 

.038 

.041 

.045 

.048 

.051 

.049 

.055 

n 

-.004 

0 

-.007 

0 

-.008 

0 

-.010 

0 

-.008 

o 

-.004 

0 

.001 

0 

.008 

0 

~76TF 

0 

~7o2T 

0 

K = — F(a/R1 ,0) 

✓5a  1 
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Upper  values  - F 
Lower  values  “ Fjj 


TABLE  10 


RANGE  OF  PARAMETER  STUDY  FOR 
FASTENER  ROW  PROGRAMS 


Fastener  hole  diameter  = 0.625  inch 
Fastener  hole  radius,  R = 0.3125  inch 


Crack  Size,  a (in.) 

s=a/ (R+a) 

a/R 

0.005 

0.0158 

0.016 

0.03 

0.0875 

0.096 

0.05 

0.138 

0.16 

0.10 

0.242 

0.32 

0.20 

0.390 

0.64 

0.30 

0.490 

0.96 

0.35 

0.528 

1.12 

94 


TABLE  11 


SINGLE  FASTENER  ROW  (CENTER  HOLE 
DAMAGED,  C/D=4) 


\a/R 

.016 

.096 

.16 

.32 

.64 

.96 

1.12 

.0158 

.0875 

.138 

.242 

.390 

.490 

.528 

90 

-.019 

0 

.020 

0 

.054 

0 

.149 

0 

.223 

0 

.304 

0 

■1 

75 

-.041 

.098 

■ 

m 

.161 

.014 

.234 

.042 

.307 

.056 

.321 

.072 

60 

-.073 

.184 

.023 

.050 

.058 

.085 

.183 

.056 

.276 

.032 

. . . 

.358 

.152 

.373 

.201 

45 

-.080 

.226 

.042 

.029 

.081 

.021 

.215 

.008 

.341 

.089 

.438 

.255 

.473 

.295 

30 

-.060 

.187 

.082 

.052 

.141 

.007 

.286 

.130 

.434 

.207 

.549 

.325 

.595 

.364 

15 

-.022 

.085 

.119 

.015 

.194 

.025 

.344 

.112 

.490 

.142 

.590 

.189 

m 

H 

m 

.149 

.006 

.228 

.010 

.367 

.012 

.521 

.038 

.611 

.048 

.656 

.054 

-15 

.088 

.013 

.168 

.042 

.242 

.090 

.367 

.151 

.546 

.254 

.625 

.298 

m 

-30 

.129 

.084 

.170 

.056 

.227 

.116 

.326 

.201 

.523 

.406 

.600 

.448 

.649 

.481 

-45 

.122 

.157 

.135 

.050 

.172 

.090 

.238 

.126 

.403 

.429 

.462 

.415 

.506 

.430 

-60 

m 

.069 

.030 

HI 

.119 

.043 

HI 

.257 

.272 

.278 

.287 

-75 

.027 

.104 

.012 

.008 

.018 

.016 

.026 

.021 

.086 

.108 

.113 

.125 

.121 
. 134 

-90 

.007 

0 

-.009 

0 

-.009 

0 

-.012 

0 

.041 

0 

.060 

0 

.070 

0 

K = — F (s , 0 ) 
/va 


Upper  values  - F^. 
Lower  values  - F 
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TABLE  12 


SINGLE  FASTENER  ROW  (CENTER  HOLE 
DAMAGED , C/D  VARIED) 


TABLE  14 


SINGLE  FASTENER  ROW  (LEFT  HOLE 
DAMAGED,  C/D=4) 


\VR 

.016 

. 096 

.16 

. 32 

. 64 

.96 

1.12 

o\S 

e X 

.0158 

. 0875 

. 138 

. 242 

. 390 

.490 

. 528 

0 

.029 

. 152 

. 232 

. 374 

. 532 

. 626 

.672 

.007 

.007 

.013 

.024 

.067 

.094 

.108 

15 

-.014 

.137 

.217 

. 378 

.538 

.648 

. 690 

.085 

.014 

. 026 

.113 

.129 

.159 

.178 

30 

-.047 

.108 

. 174 

. 337 

. 508 

.633 

. 685 

.185 

. 046 

.016 

.155 

. 230 

. 321 

. 351 

.063 

.108 

.255 

.406 

. 516 

. 561 

1 

.034 

.046 

.160 

.286 

. 315 

60 

. 030 

BSi 

.199 

. 303 

.404 

.429 

.178 

.060 

HU 

.030 

.033 

.191 

.232 

75 

-.040 

. 019 

. 056 

.154 

.236 

. 324 

. 350 

.096 

.093 

.122 

.004 

.001 

.103 

. 120 

90 

-.020 

.016 

.143 

.222 

.317 

.353 

0 

0 

.012 

.030 

.044 

.049 

105 

-.041 

.017 

.052 

.156 

.239 

. 322 

. 342 

.099 

.086 

.122 

.027 

.080 

.013 

.024 

120 

-.078 

.012 

.044 

. 168 

.269 

. 354 

. 366 

.188 

.039 

.084 

.104 

.115 

.170 

135 

mam 

.290 

. 382 

.412 

8H 

.012 

.235 

.286 

150 

-.072 

. 052 

. 104 

. 234 

.368 

. 524 

.188 

.061 

.001 

.105 

.193 

m M 

.416 

165 

-.029 

.101 

.172 

. 312 

.451 

.555 

. 601 

.086 

.015 

.027 

.116 

.170 

. 252 

.298 

K = — F (s,  0) 
/Sa 


Upper  values  - F^ 
Lower  values  - F^ 
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TABLE  14 


(Concluded) 


\*/R 

\« 

.016 

.096 

.16 

.32 

.64 

.96 

1.12 

9 \ 

.0158 

.0875 

.138 

.242 

.390 

.490 

.528 

180 

.029 

.153 

.234 

.375 

.535 

.634 

.685 

.007 

.006 

.007 

. 001 

.011 

. 009 

.009 

195 

.095 

.189 

. 268 

.404 

.597 

. 685 

.725 

.009 

.046 

.095 

.163 

. 268 

. 304 

. 336 

210 

.137 

.190 

.252 

.361 

.568 

.639 

.685 

.077 

.069  , 

.133 

.247 

.480 

.503 

.534 

225 

.126 

.141 

.179 

.242 

.402 

.446 

.485 

.152 

.074 

.068 

.111 

.183 

.541 

.490 

.500 

240  | 

.062 

. 0 79  1 

.098 

.170 

.201 

.215 

.163  ! 

f 

.042 

.064 

.079 

i 

.348 

- ■ - 

.324 

.336 

! 255  ! 

L J 

.025  | 

.001  | 

.002 

-.001  i 

.040  I 

i • 054 ; 

.056 

.100  j 

.013  i 

l 

.025 

.039  j .144 

. :162i -173_ 

270  | 

) 

l 

- 006  T 

-.018  | 

-.023 

-.036  | 

.002 

.011 

.015 

0 

.002 

.007  j 

.014 

.018|  .019 

285  ! 

j 

.023  ! 

.003 

.004  ! 

.004  i 

.052 

.070 

.075 

.101  i 

| 

.008  | 

. 018  j 

.020  i 

i 

.099 

.121 

.133 

300  j 

.068  ; 

.055  j 

.069  j 

.094  j 

.174 

.210 

. 224"” 

.163  ! 

J 

.025  | 

.047s 

\ 

.034  , 

.254 
. 349’ 

.277 

.300 

j 

315  j 

i 

.no  i 

.113  j 

.144 : 

.199 

.402!  .439  1 

.155  ! 

.042  i 
' .146'! 

. 085  j 

.100  i 

. 393 

.423 

.448 

r 

330  : 

| 

. iie] 

. 199  j 

. 286  I 

. 470 

.541 

. 590 

. 085  j 

.051 

.113! 

.185  | 

.405 

.477 

.520 

1 

34  5 j 

.081  j 

.155 

.226; 

. 345  j 

.516 

.595 

.628 

.015  ! 

.042  j 

.092! 

.155  1 

i 

.277 

. 342 

.351 

K = — F (s, 0) 
/ira 


Upper  values  - F 
Lower  values  - F 
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TABLE  15 


DOUBLE  FASTENER  ROW  (CENTER  HOLE 
IN  UPPER  ROW  DAMAGED,  C/D=4) 


\*/R 

\ 

\ 

.016 

.096 

1 

i 

.16 

.32 

. 64 

l I 

.96  i 

1.12  j 

\ c 

a”  Xs 
\ 

| 

.0158 

i 

.0875!  .138 

. 242: 

. 390 

.490; 

. 528  1 

0 

.029 

.011 

.155 

.006 

.237 
.0  07 

^ .382! 

.003; 

.543 

.014 

. 632! 
.009  ' 

.675  ! 
.010  ; 

— 

15 

pToi51 

.110 

.119 

.013 

.195 
.0  28 

. 354 
.125' 

1 

j _ 

.490 

.158 

. 589 
, .218, 

.626  ! 
.251 

30 

nr;  103 

. 259 

| .081 
.053 

.141 

.006 

.302: 
.14  0: 

i - 

.432 

.186 

.554; 

.311: 

_ 

.602  ; 
.349  j 

45 

I-.122 
. 338 

.045 

.033 

.085 

.015 

. 

.241* 

.017; 

. 362 
.034 

.479 
. 208; 

.522 
.242  ! 

60 

'-.102 

.303 

.028 

.046 

.067 

.077 

.216! 

.048 

. 315 
.085 

.427, 

.097: 

.456~1 
.136  : 

75 

-.055 
1 .176 

.027 

.086 

.069 

.109 

.194! 
; .012 

. 275 
.073 

1 .380; 

.018; 

.410 
.025  ■ 

90 

-.0  25 
0 

.030 

^7070 

.001 

! .188! 
j .004; 

. 271 
.009 

j .386j 
| .013| 

.4  29  : 
.015  j 

105 

! - . 0 50 

i .175 

1 ! 

.038 

.085 

! .082 
| .109 

■ .213: 

: .012; 

j - 

.301 

.065 

; .409' 

, . 038, 

.439  ! 
.050  . 

120 

-.095  i 
i .303  | 

.044 

.041 

' .086 
.072 

: .242: 

; .064! 

1 j 

.351 

.102 

.466: 

.098 

. 495  ; 
.144  j 

135 

— .116  ■ 
I .340  : 

.057 

.040 

i .loo 
; .oo7 

. 261 
.007' 

.387 

.012 

i .505: 
.189; 

.550  : 
.229  ; 

i 

150 

1 

-.100  1 
.260 

.086  j 

.059 

.147 

.001 

.309* 
.12  5' 

1 T 

.437 

.151 

: . 561} 

! .292; 

.611 
.333  ; 

1 

i 

165 

-.050 

.110 

.120  ! 
.015 

.196  ! 
.026 

.355! 
. 118: 

.487 

.138! 

! .588! 

1 •1971 

.625  j 
.229  j 

p 

K - F(s,0)  Upper  values  - 

Lower  values  - F „ 
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TABLE  15 


(Concluded) 


\a/R 

.016 

.096 

.16 

.32 

.64 

— 

• | 

VO  ! 

cn  ; 

. J 

1.12 

o\  c 

9 

.0158 

.0875 

.138 

.242 

.390 

.490 

.528 

180 

.029 
. Oil 

.155 

.006 

.237 

.009 

.38  2 
.004 

.542 

.029 

.632 

.032 

.674 

.034 

195 

mn 

■ 

.180 

.043 

.259 

.092 

. 385 
.148 

.586 

.255 

.666 

.284 

.700 

.311 

210 

.118 

.153 

.185 

.060 

.247 

.124 

.343 

.203 

.573 

.430 

.652 

.462 

.708 

.492 

225 

m 

.146 

.059 

mm 

.245 

.128 

.449 

.498 

.514 

.465 

.565 

.480 

240 

.108 

.286 

.071 

.038 

mm 

ESI 

• 111 
.045 

.229 

.346 

.279 

.330 

.304 

.351 

255 

m 

.005 

.012 

.009 

.023 

.004 

.023 

.079 

.146 

.099 

.164 

.104 

.180 

270 

.014 

0 

-.015 

0 

-.016 

.001 

-.032 

.004 

.032 

.009 

.044 

.014 

.052 

.017 

285 

.044 

.181 

.016 

.011 

.023 

.023 

.024 

.022 

■Bill 

1 ' ' 

.128 

.185 

.135 

.207 

300 

.114 

.288 

.086 

.034 

.110 

.056 

.139 

.029 

.267 

.330 

.320 

.333 

.348 

.360 

315 

.176 

.270 

.159 

.051 

.202 

.095 

.266 

.104 

.476 

.452 

.575 

.446 

.598 

.467 

330 

HB9 

■SI 

.190 

.055 

.254 

.118 

.351 

.188 

.578 

.393 

.660 
. 443 

.718 

.476 

345 

.120 

.036 

.181 

.042 

.260 

.090 

.387 

.141 

.585 

.235 

.665 

.264 

.700 

.289 

K = — F(s,9) 

/ ia 


Upper  values  - F 
Lower  values  - F 
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TABLE  16 


DOUBLE  FASTENER 
IN  LOWER  ROW 


ROW  (CENTER  HOLE 
DAMAGED,  C/D=4) 


\a/R 

.016 

.096 

.16 

.32 

.64 

.96 

1.12 

0 \s 
e \ 

■H 

■Aififiii 

.0875 

.138 

. 242 

.390 

.490 

. 528 

0 

.046 

.009 

. 243 
.007 

.374 

.012 

.604 

.019 

. 872 
. 062 

1.036 

.082 

1.119 

.093 

15 

. 342 
. 076 

.586 

.198 

.841 

.265 

1.018 

.355 

1.092 

.410 

30 

-.044 

.261 

.187 

.024 

.284 

.080 

. 514 
.258 

.763 

.437 

.954 

.646 

1.040 

.722 

45 

-.066 

.317 

.128 

0 

. 192 
.088 

. 390 
. 120 

.589 

.377 

.744 

.616 

.811 
. 682 

60 

-.071 

.275 

.061 

.076 

.263 

.036 

.373 

.214 

.494 

.464 

. 519 
.545 

El 

.009 

.101 

. 039 
.145 

m 

.208 

.087 

.294 

.238 

90 

-.033 

0 

-.016 

0 

. 003 
.001 

m 

.138 
. Oil 

. 220 
.016 

105 

-.057 

.160 

WSmWivm 

■fiMfl 

.022 

.145 

.123 

.054 

.176 

.079 

.259 

.215 

.270 

.247 

120 

-.083 

.270 

.043 

.081 

.083 
. 144 

.232 

.056 

.331 

.241 

.450 

.466 

.474 
. 589 

135 

-.076 
. 307 

.115 

.010 

.175 

.098 

. 371 
.150 

. 563 
.438 

.784 

.706 

150 

-.049 

.250 

.182 

.017 

.278 

.087 

Baa 

.763 

.484 

. 956 
. 675 

1.042 

.749 

165 

bob 

.220 

.009 

. 342 
.080 

. 588 
.206 

. 849 
. 291 

1.030 
. 385 

1.107 

.442 

K = — F (s , 6 ) 
/rra 


Upper  values  - 
Lower  values  - F^ 
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TABLE  16 


(Concluded) 


TABLE  17 


DOUBLE  FASTENER  ROW  (CENTER  HOLE 
IN  LOWER  ROW  DAMAGED  AND  UNLOADED,  C/D=4) 


o \s 

r -N 

0 

.016 

.096 

1 .16 

.32 

. 64 

.96 

1.12 

I 

.0158,  .0875 

. 138 

. 242 

. 390 

. 490 

. 528 

. 024 
. 006 

.134 

0 

.209 
. 001 

. 348 
0 

. 527 
.019 

. 654 
. 021 

. 721 
.023 

-.009 

.140 

. 209 

. 354 

. 517 

. 649 

. 708 

15 

. 060 

.030 

. 072^ 

.125 

. 179 

. 250 

.290 

-.015 

.149 

.203 

. 339 

.486 

. 619 

.687 

30 

. 142 

.037 

. 098 

. 181 

. 309 

. 453 

. 512 

-.017 

. 127 

.162 

.274 

. 379 

. 481 

. 534 

4 5 

.200 

.036 

. 086 

.157 

. 365 

.485 

. 529 

- .02  3 " 

.063 

. 078  ' 

.146 

.175 

. 245 

. 266 

6 0 

.208 

. 030 

.063 

.127 

. 300 

. 401 

. 446 

-.024 

-.005 

-.010 

. 014 

. 005 

.041 

. 045 

! / 5 

> 

.144 

. 014 

.034 

. 091 

.153 

.227 

. 256 

n r\  ; 

-.023  ’ 

-.040  ! 

-.056 

- . 0 57 

-.082 

061 

-.059” 

y 0 ! 

i 

j 

0 

0 

. 001  j 

. 004 

. 012  j 

.018 

.021 

1 1 

105  1 

1 

1 

- . 027  ■ 

-.026  : 

-.037  : 

-.026  1 

-.049  | 

- . 019' 

- . 0l.6~ 

.150  ; 

. 018  , 

. 037  : 

.100  j 

.157 

. 206 

. 225 

1 

j 

-.026  ‘ 

. 0 39  i 

. 045  j 

. 095  j 

. 108  j 

* .176” 

.195" 

120  | 

1 

. 224  j 

. 045  ; 

i 

. 080  j 

.179  j 

. 374  ! 

.422 

.452 

1 

1 

-.014  ! 

. 115  : 

.147  | 

. 253 

.358  j 

.4  6 0 ~ 

."512” 

135  ! 
1 

.221  j 

. 060  1 

. 110  | 

. 239 

. 508 

. 551 

. 579 

j 

i r a { 

-. 008  ! 

.155  ] 

.210  ; 

’”.353" 

.519  " 

.6  55' 

".  723” 

150  | 

. 164  i 

i 

.053 

.114  ; 

. 234 

. 401 

. 501 

. 549 

t rr 

-.003  | 

. 150  J 

.221  1 

. 372 

. 5 53 

. 691 

.754 

165 

.076 

.034  j 

. 078  | 

. 145 

. 214  j 

. 280 

. 318 

p 

K = — F(s,0)  Upper  values  - F 

/«T  ra 

‘ c'  Lower  values  - F^ 
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TABLE  17 


(Concluded) 


V/k 

.016 

.096 

.16 

.32 

.64 

Xf 

e X 

.0158 

.0875 

.138 

.242 

.390 

.490 

.528 

180 

■ 

.136 

0 

.212 

0 

.351 

.006 

.539 

.004 

.672 

0 

.742 

0 

195 

B 

m 

.207 

.074 

.326 

.113 

.531 

.202 

.653 

.272 

.712 
. 317 

210 

.131 

.078 

.146 

.037 

m 

.291 

.140 

. 511 
.333 

225 

.130 

.156 

.125 

.034 

.160 

.082 

.223 

.079 

i 

.407 

.384 

.486 

.462 

.529 

.503 

240 

.085 

.174 

.064 

.028 

.079 

.058 

m 

.206 

.316 

.257 

.368 

warn 

■iliKl 

Wmmmm 

255 

HSQ| 

BHWJH 

-.005 

.030 

-.012 

.043 

.039 

.155 

.064 

.197 

.071 

.219 

270 

EEBI 

BBM 

■□i 

-.047 

.001 

-.074 

.005 

-.040 

.012 

m 

m 

285 

EES 

1 

a 

m 

-.009 

.151 

.011 

.171 

m 

300 

.064 

.187 

.037 

.040 

. __ 

.046 

.070 

.059 

.082 

.135 

.365 

.182  1 
.369 

.201 

.393 

315 

.113 

.181 

m 

.140 

.100 

■g 

mm 

.361 

.495 

.434 

.501 

.481 

.526 

330 

IBB 

.145 

.051 

■BEE1 

1 

.293 

.174 

.508 

.411 

.614 

.484 

.674 

.530 
! 

345 

m 

.143 

.034 

• 212 
.078 

.335 

.124 

.538 

.234 

.655 

.294 

.711 

.336 

Upper  values  - F^ 
Lower  values  - 


K = — F(s,6) 

/FI 
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